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Abstract. The long-term variability of a sample of 180 
optically selected QSOs in the field of the Selected Area 
94 has been studied. The relations between variability and 
luminosity and between variability and redshift have been 
investigated by means of statistical estimators that are 
"robust" and allow at the same time to eliminate the in- 
ON fluence of the measurement errors. A comparison is carried 
out with the results of two other samples of QSOs, in the 
fields of the South Galactic Pole (Hook et al. 1994) and 
of the Selected Area 57 (Trevese et al. 1994). Merging the 
three samples provides a total of 486 QSOs. The analysis 
in the QSOs rest frame of both the ensemble structure 
function (SF) and the individual variability indices show 
that: 1) a negative correlation between variability and lu- 
minosity is clearly present, in the sense that more lumi- 
nous QSOs show less variability; 2) a significant positive 
correlation exists between variability and redshift; 3) such 
correlations may be equally well parameterized either with 
a model in which the timescale of the variability is fixed 
for all the QSOs (r ~ 2.4 yr), while the amplitude lin- 
"-^ early increases with the absolute magnitude and redshift, 
or with a model in which the timescale of the variabil- 
ity linearly depends on the absolute magnitude and the 
amplitude is only a function of the redshift. 

The same analysis carried out in the observer's frame 
^ provides the following results: 1) there is a negative corre- 
lation between variability and luminosity; 2) the timescale 
of variability does not depend significantly either on the 
absolute magnitude or on the redshift; 3) the ensemble 
structure function is well represented by a parameteriza- 
tion in which, with a fixed timescale of about 5.5 yr, the 
amplitude linearly increases with the absolute magnitude; 
4) although the general behaviour of the SF does not show 
a systematic variation of the timescale and/or amplitude 
with redshift, if we examine the average variability index 
for objects with —25 > Mb > —27, we find that below 
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* Based on material collected with the UKSTU and the ESO- 
La Silla telescopes and on COSMOS scans. 



redshift 1 quasars are significantly less variable than at 
higher redshift. 

The implications in terms of the black-hole, starburst 
and microlensing models are briefly discussed. 
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1. Introduction 

Variability is a common property of QSOs. Although 
not thoroughly understood and exploited, it provides a 
potentially powerful tool to constrain theoretical mod- 
els. A number of different scenarios, such as the black- 
hole (Rees 1984), the starburst (Terlevich et al. 1992), 
the gravitational microlensing (Hawkins 1993) models, 
in principle can be tested on the basis of their predic- 
tions about the temporal variations of the QSO lumi- 
nosity. Attempts in this sense have been reviewed by 
Wallinder et al. (1992) in the framework of the black- 
hole model. Up to now, unfortunately, the results have 
not come up to general expectations. Even the determi- 
nation of simple phenomenological trends, such as the 
correlations of variability with redshift or luminosity, 
has resulted in a complex situation (Cristiani et al. 1990, 
Giallongo et al. 1991, Hook et al. 1994, Trevese et al. 
1994, Paltani &; Courvoisier 1994). 

Other important results from the study of the ensemble 
variability of QSOs concern the understanding of the bi- 
ases affecting samples selected on the basis of the sole vari- 
ability (Trevese et al. 1989, Veron &: Hawkins 1995) or of 
multicolor data (Warren et al. 1991). The analysis of the 
relationship among the luminosity in different bands of 
the QSO e.m. spectrum is also affected by variability (e.g. 
the L x — L a relationship, La Franca et al. 1995). 

Variability provides an easy way to select independent 
QSO samples to cross-check the biases of various selection 
techniques (Cimatti et al. 1993, Trevese et al. 1994). 



standing at least of the above mentioned phenomenolog- 
ical trends prompted us to continue and refine the study 
of the field of the Selected Area 94 (Cristiani et al. 1990), 
developing statistical techniques allowing a meaning- 
ful comparison with the other major variability studies 
(Hook et al. 1994, Trevese et al. 1994). 

In the following we report in Section 2 the description 
of the QSO of the Selected Area 94 (SA 94) sample and 
the photometric material used to study its variability; in 
Section 3 the characteristics of two other samples in the 
Selected Area 57 (SA 57) and South Galactic Pole (SGP) 
are recalled; in Section 4 the ensemble structure function 
is introduced as a statistical measure of the variability 
properties and the detected trends are further analysed 
with the aid of another variability index; the results and 
implications are discussed in Section 5. 

Throughout the paper we have assumed cosmological 
constants H a = 50 Km sec -1 Mpc -1 , q a = 0.5. 

2. The SA94 sample 

2.1. Definition of the sample 

The SA94 sample is made up of 180 optically selected 
QSOs 1 (listed in Table 1), observed with 21 plates cover- 
ing a time-base of 10 years. The QSOs belong to the UVx 
and objective-prism surveys carried out by Cristiani et al. 
(1991), La Franca et al. (1992), Cristiani et al. (1995); who 
also discussed the properties of the surveys in terms of se- 
lection criteria, completeness (better than 90%) and detec- 
tion efficiency. 47 previously unpublished QSOs have been 
included in Table 1 in order to enlarge the sample and in- 
crease the statistical significancy. They have been selected 
with the same criteria from the same databases as the rest 
of the sample and therefore are expected to share the same 
variability properties. Their identification has been carried 
out with spectroscopic observations obtained at the ESO 
3.6-m telescope in La Silla, equipped with the OPTO- 
PUS multifibre spectrograph (Avila &: D'Odorico 1993), 
in three observing runs on September 1990 and Septem- 
ber and October 1991. The OPTOPUS multifiber spectro- 
graph uses bundles of 50 optical fibers, which can be set 
within the field of the Cassegrain focal plane of the tele- 
scope; this field has a diameter of 32 arcmin, and each fibre 
has a projected size on the sky of 2.5 arcsec. We used the 
ESO grating #13, providing a dispersion of 450 A mm -1 
in the wavelength range from 3800 to 9800 A. The res- 
olution was 27 A. We dedicated about half of the fibers 

1 In the present paper QSOs are defined as objects with a 
starlike nucleus, broad emission lines, brighter than Mb = — 23 
mag, applying the K-corrections computed on the basis of the 
composite spectrum of Cristiani and Vio (1990), and dered- 
dened for galactic extinction according to Burstein and Heiles 
(1982). 



Table 1. List of the QSOs in the SA94 
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to the sky. The observing time for each field was 1 hour, 
split into two half-hour exposures performing object-sky 
flipping. A S/N ratio per resolution element larger than 10 
was obtained, allowing clear identification for more than 
90 per cent of the QSO candidates. 

2.2. Calibration of the photographic material and error es- 
timation 

21 plates taken with the ESO La Silla and UK Schmidt 
telescopes have been analysed. The combination Ila-O + 
GG385 filter was always used, defining a passband close to 
the Johnson B. The color transformation between these 
"natural" B' photographic magnitudes and the Johnson 
B is (Blair & Gilmore 1982) 

B' = B - 0.11(5 - V). (1) 

The field of the sky investigated is included into the 
limits of the Selected Area 94 (1950.0 coordinates in the 
following intervals: 2 ft 43™27.5 s < a < 3 ft l™34.0 s and 
-2°05'8.9" < 8 < 2°46'51.2"), covering an area of 22.03 
square degrees. In Table 2 a detailed list of plates is given 
(the meaning of limit- mag is defined below). 

The plate material has been scanned with the COS- 
MOS microdensitometer (MacGillivray &: Stobie 1984). 

The resulting tables, one per each plate, containing the 
instrumental magnitudes and other useful parameters for 
the objects detected, have been merged together in one 
table. Only objects with at least 4 detections in the first 
10 plates have been accepted in this final table. The as- 
trometric error box individuating a common detection has 
been defined as a circle of 3.5 arcsec of radius. In this way 
spurious detections (plate defects) are minimized to an 
acceptable level, while real measurements are in practice 
never discarded. 

To put all the plates in a common magnitude scale: 

1) the instrumental magnitudes of each plate sepa- 
rately have been calibrated using polynomial regressions 



Table 2. List of the Plates 



Plate Exp Date Epoch/date Limit 



(min) (mag) 



B7221 


60 


1981 


Sep 28 


1 


(1981 


74) 


20 


3 


B5411 


60 


1983 


Dec 07 


2 


1983 


94) 


21 


1 


B5415 


60 


1983 


Dec 08 


2 


1983 


94) 


20 


9 


B6321 


60 


1986 


Jan 05 


3 


1986 


02) 


20 


9 


B6325 


60 


1986 


Jan 06 


3 


1986 


02) 


20 


9 


B6673 


60 


1986 


Oct 25 


4 


1986 


82) 


21 


1 


B6692 


60 


1986 


Nov 21 


5 


1986 


89) 


21 


1 


B6983 


60 


1987 


Aug 20 


6 


1987 


65) 


20 


9 


B6984 


60 


1987 


Aug 20 


6 


1987 


65) 


20 


9 


B6989 


60 


1987 


Aug 24 


6 


1987 


65) 


20 


9 


B6990 


60 


1987 


Aug 24 


6 


1987 


65) 


20 


9 


B6993 


60 


1987 


Aug 25 


6 


1987 


65) 


20 


5 


B6994 


60 


1987 


Aug 25 


6 


1987 


65) 


20 


5 


B7177 


60 


1987 


Dec 20 


7 


1987 


97) 


20 


9 


B7703 


60 


1988 


Nov 28 


8 


1988 


91) 


20 


9 


B8268 


75 


1989 


Oct 24 


9 


1989 


82) 


20 


7 


B8280 


75 


1989 


Oct 27 


9 


1989 


82) 


20 


9 


B8742 


70 


1990 


Sep 14 


10 


(1990 


.71) 


20 


7 


B8746 


70 


1990 


Sep 21 


10 


(1990 


.71) 


19 


9 


B9405 


75 


1991 


Aug 13 


11 


(1991 


.62) 


20 


7 


B9523 


75 


1991 


Oct 31 


12 


(1991 


.83) 


20 


7 



obtained via 116 photometric standards available in the 
literature or observed specifically for this programme, cov- 
ering the interval 6.5<5<21.2. 

2) for each object for which at least 5 measurements 
were available, the median of the calibrated magnitudes 
has been computed, defining a set of reference magnitudes 
{B' ref , 48 863 objects). 

3) to obtain a useful parameter for the separation of 
extended from point-like sources, we have computed for 
each plate and each object the difference between the mea- 
sured FWHM and the mode of the FWHM distribution at 
the magnitude of the object. This quantity has then been 
averaged for each object on the first ten plates. 

4) for each plate, the magnitudes computed in step 1) 
have been re-calibrated against the reference magnitudes, 
using only point-like objects (41 151 objects). 

5) a procedure of uniformization of the usually spa- 
tially variable response of the photographic plates has 
been applied. Each plate has been subdivided in 10 x 10 
sub-areas, for each of them the differences between the 
reference and the individual plate magnitudes have been 
computed and their distribution analysed. The zero-point 
shifts estimated in this way for each sub-area as a function 
of the magnitude have been smoothed and applied to the 
re-calibrated magnitudes. In the following we will refer to 
the magnitudes obtained in this way as B'^ inal . 

For each plate the limit of completeness has been as- 
sumed to coincide with the maximum of the histogram of 



last column of Table 2). 

For each plate the uncertainties on the B'^ inal magni- 
tudes have been estimated by analysing the distribution 



of the differences AB = B' t 



final 



B' re f as a function of the 



reference magnitudes for all the point-like objects. 

3. The SA57 and SGP samples 

The SA57 sample is made up of 23 QSOs, studied by 
Trevese et al. (1994) using 14 plates covering a time-base 
of 15 years. To evaluate the magnitude errors pertaining 
to each object we have applied a regression to the data of 
Table 4 of Trevese et al.'s paper, showing the dependence 
of the photometric error as a function of the apparent 
magnitude. For the structure function and variability in- 
dex calculation (see below) we have considered the same 
subdivision in 11 epochs applied by Trevese and collabo- 
rators. 

The South Galactic Pole sample, comprising 283 
QSOs, has been studied by Hook et al. (1994) using 11 
plates, covering a time interval of 16 years. We have subdi- 
vided the 11 plates in 7 epochs, the same as in the original 
paper (Hook et al. 1994). To evaluate the magnitude er- 
rors, the same procedure as in the case of the SA57 sample 
has been applied (see Table 2 of Hook et al.'s paper). 

4. Statistical indices of the variability 

4-1. The structure function 

One of the simplest and most immediate methods to 
analyse the variability of an object is to calculate the 
structure function (SF). Properties and limitations of 
the SF and its application to derive timescales and am- 
plitudes of the QSO variability have been described by 
Simonetti et al. 1985, Trevese et al. 1994, Vio et al. 1992 
and references therein. For a finite sequence of measure- 



N: 



ments f(t), t = 1,2. 

SF{T) = Np)^< W{t)w{t + T)[m{t + T) ~ m{t)]2 (2) 

where m(t) are the magnitudes at the time t, w(t) are the 
weights, equal to 1 if a measure exist for the time t, 
otherwise, and 



N(t) = ^2w(t)w(t + r) 



(3) 



In the present case it is not possible to study in detail 
the variability of the individual QSOs, due to the relatively 
small number of plates. We can instead treat the objects 
in our sample as representative of the QSO "class" and 
compute a quantity analogous to the SF: 



SF e {r) = ([mag{t) - mag{t + r)] 2 ) 



(4) 



where the r is evaluated in the rest-frame of each QSO and 
the brackets " (. . .) " indicate a mean over the ensemble. 



of this approach. 

Equation (4) does not provide an optimal estimation 
of the SF in the sense of statistical "robustness" . It de- 
pends on the second moment of the distribution of the 
magnitude differences, which is exposed to the influence 
of possible "outliers" . In principle an estimator based on 
the mean absolute deviation is expected to give better re- 
sults, but the presence of the measurement errors (whose 
subtraction is a necessary condition to obtain a "sample- 
independent" SF) makes its computation complex. To 
this end, we adopted the procedure described in the ap- 
pendix. 
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Fig. 1. Upper panel: structure functions for the QSO samples 
SA94 (filled circles), SGP (crosses) and SA57 (open squares). 
The open circles show the structure function evaluated on the 
control sample (see text). Lower panel: global structure func- 
tion obtained from the three QSO samples 



Different scenarios require different types of computa- 
tion of the SF. If the variability is assumed to be intrinsic, 
as in the black-hole or the starburst model, it has to be 
evaluated in the QSOs restframe: time intervals between 
pairs of epochs have to be scaled by the (1 + .z ?so ) _1 cos- 
mological factor. On the other hand, Hawkins (1993) has 
recently suggested an "intervening" origin for the QSO 
variability, with nearly all QSOs being microlensed. Ac- 
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Fig. 2. QSO distribution in the M B - z plane: SA94 (filled 
circles), SGP (crosses), SA57 (open squares). Absolute magni- 
tudes have been calculated according to Cristiani fc Vio (1990) 
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cording to him the redshift distribution of the lensing ob- 
jects (for a uniform population of lenses) is expected to 
strongly peak at z ~ 0.5, with only weak dependence on 
the redshift of the source, and therefore the observer's 
rest-frame is the natural choice for the SF computation. 



Fig. 3. The structure function for QSOs of different luminosity 
in the redshift interval 0.3 < z < 1.4. Upper panel: QSOs with 
M B < -24.5. Middle panel: QSOs with -24.5 < M B - Lower 
panel: comparison between the two SFs. Symbols are the same 
as in Fig.l 



4.1.1. The structure function in the QSOs rest frame 

In the lower panel of Fig. 1 the rest frame SF computed 
simultaneously over the SA94, SGP and SA57 samples, 
with their error-bars, is shown. We observe a steady rise, 
with a tendency towards flattening, reaching values cor- 
responding to about 0.1 mag 2 . We stress here that the 
shape of the ensemble structure function is the result of 
the combination of the intrinsic variability properties of 
the QSOs that contribute in different proportions to each 
bin of time (for example the larger At's are populated 
mainly by lower-redshift objects) and is truly representa- 
tive only if all the QSOs in the samples behave more or 
less in the same way. This is not the case, as will be shown 
below. 

In the upper panel of Fig. 1 the SFs for the SA94, 
the SA57 and the SGP samples are shown. The structure 
function for a control sample in SA94, made of known stars 
and point-like objects lacking a spectroscopic classification 
(41 151 objects), is reported too. While the SFs for the 
SA94 and SGP samples are very similar for t < 3yr, they 
diverge for longer periods of times. Besides, the QSOs in 
SA57 present a remarkably larger SF for timescales t < 



3yr. As indicated by previous works (Cristiani et al. 1990, 
Hook et al. 1994, Trevese et al. 1994), differences in the 
SFs may be due to correlations between the variability 
and physical parameters like absolute luminosity and/or 
redshift, combined with a different coverage of the L — z 
plane in the various samples. In Fig. 2, the distribution 
of all QSOs in the redshift-magnitude plane is illustrated. 
We see that SA57 objects are, on average, less luminous 
than the others. 

To further investigate and disentangle the above men- 
tioned dependences we have examined two subsamples. 
The first one is defined within the redshift limits 0.3 < 
z < 1.4 (the area between the dashed lines in Fig. 2). The 
SF e of less luminous QSOs {Mb > —24.5) shows a larger 
amplitude (Fig. 3) with respect to the more luminous ones 
{Mb < —24.5). The second subsample is defined within 
the absolute luminosity limits —26.85 < Mb < —25.5 (the 
area between the continuous lines in Fig. 2). The SF e for 
low-redshift QSOs {z < 2) is not distinguishable, within 
the errors, from the one of high-redshift QSOs {z > 2), as 
shown in Fig. 4. No obvious change is apparent neither of 
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the amplitude nor of the timescales of the variability as a 
function of the redshift. 

To further investigate the dependences of the SF on 
the absolute magnitude and redshift we have subdivided 
the L — z plane in a number of sub-areas and evaluated 
the SF e in each of them. The result is shown in Fig. 5. 

Going from upper to lower panels within the various 
redshift slices, the anti-correlation between the SF ampli- 
tude and absolute magnitude is apparent. A much weaker 
dependence (if any) exists on the redshift. 

In order to parameterize the observed properties of 
QSO variability with luminosity, redshift and rest-frame 
time interval, we considered some simple functional forms 
in which the dependences are separable. Two basic forms 
of the SF have been explored, a negative exponential 
SF = A(l - e-'/ T ) and a power-law SF = At T . The 
amplitude A and the timescale r have been in turn as- 
sumed to be a linear function of the absolute magnitude 
or redshift, producing six different models. 



SF = [a + b(M B + 25.7) + cz]{\ - e~ t/T ) (5) 
Model B: 



SF = [a + b(M B + 25.7) + cz]t T 
Model C: 

SF = [a + b{M B + 25.7)](1 - e -*/^+ c (i+^)]) 
Model D: 

SF = [a + b{M B + 25.7)]^ T+C ( 1+2 )] 
Model E: 

SF = (a + CZ)(1 - e -t/[r+6(M B +25.7)]j 

Model F: 



SF = (a - 



■ cz)i [T+h( - MB+2Z - 7)] 



(6) 



(7) 



(8) 



(9) 



(10) 



Each model was used to fit the binned data shown in 
Fig. 5 with a non-linear least-squares method. The bin 
at < M B >= -24.3 and < z >= 2.28, based on only 4 
objects pertaining to the SA57, appears highly discrepant 
from the general trend and has been excluded from the 
fitting procedure (see the next subsection for further com- 
ments on this point). The results are given in Table 3. 

Errors quoted for the parameters are 68 per cent con- 
fidence intervals. Model A and E give the best fit to the 
data. In both cases dropping the dependence on the red- 
shift provides a significantly worse result. The F-test on 
the inclusion of the parameter c gives respectively a 5-10 -3 
and a 2 • 10 -2 probability that the improvement obtained 
including the z-dependence is due to chance. 

4.1.2. The structure function in the observer's rest frame 

The same parameterizations of the SF have been inves- 
tigated in the observer's frame, to carry out a compari- 
son with the microlensing model of the QSO variability 
(Hawkins 1993, Lacey 1994, Alexander 1995). The results 
are given in Table 4. 

Now model A and C give the best fit to the data. In 
both cases the dependence on the redshift is very weak and 
dropping it does not change the goodness of fit (model A 
and C coincide when the dependence on z is dropped). 
Thus, when analysing the variability in the observer's 
frame, the best fit is obtained with only three parame- 
ters. The amplitude of the SF depends on the absolute 
magnitude, while the timescale seems to be independent 
on absolute magnitude or redshift. 
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Fig. 5. The SF (continuous line) as a function of absolute magnitude and redshift. In the upper left corner of each panel the 
average absolute magnitude, the average redshift and the number of objects in the bin are shown. The dashed line corresponds 
to the best-fit model A of Table 3 (see text). The bin at < Mb >= —24.3 and < z >= 2.28 is re-plotted in the upper right 
corner with a lower scale in order to show all the points 



Table 3. Best-nt parameters for the functional forms A, B, C, Table 4. Best-nt parameters for the functional forms A, B, C, 
D, E and F described in the text, evaluated in the QSO rest D, E and F described in the text, evaluated in the observer's 
frame frame 
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photometric errors, see appendix) required to reproduce 
the average absolute differences between pairs of epoch- 
magnitudes (cfr. eq. A6). The difference with the previous 
4-2. The variability index section is that now the intrinsic variance is estimated for 

each object individually and not at one time collectively 
Although the union of the 3 samples has increased the for all the objects contributing to a given temporal bin. 
number of objects, allowing us to carry out tests on sub- Only pairs of epochs separated by more than 1 and less 
samples, the distribution of the objects in the L—z plane is than 4 years in the QSOs rest frames have been considered, 
far from uniform, leaving the possibility for spurious cor- in order to avoid on the one hand time scales for which the 
relations between variability and redshift to be induced structure function is still quickly rising, as observed in the 
by the still surviving correlation between redshift and ab- previous section, and on the other hand time lags unacces- 
solute luminosity. In order to clarify this point we have sible with the present data for the higher-redshift QSOs. 
followed a second approach based on the definition of a The choice of this time interval is in any case not criti- 
variability index for each QSO. By analogy with the SF, cal, and the following results depend very weakly on it, as 
we have computed a variability index (IDX) for each shown by extensive tests. The use of an index defined in 
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Fig. 6. Variability index versus apparent magnitude for the 
QSOs in the SA94 (•), SGP (x) and SA57 (□) 



this way has an immediate interpretation in terms of a pa- 
rameterization of the type of Model A (shown in the pre- 
vious sub-section to be a satisfactory representation of the 
data), for which the timescale of the rest-frame variability 
is the same for all the objects and the amplitude varies 
as a function of absolute magnitude and redshift. Even in 
the framework of a different parameterization (e.g. Model 
E) this variability index can still give useful indications, 
although of less immediate interpretation. 

In Figs. 6, 7, 8 the variability indices vs. apparent mag- 
nitude, absolute magnitude and redshift for the QSOs of 
SA94, SA57 and SGP are shown. The variability indices 
of the objects in the SA94 are reported in Table 1. 

The computation of the variability index for the con- 
trol sample, defined in the previous sub-section, shows 
that, as expected, no correlation is present between vari- 
ability and apparent magnitude, confirming that the mea- 
surement errors have been effectively removed. 

To quantify the visual impressions given by Figs. 6, 7 
and 8 we have analysed the correlation matrix, reported 
in Table 5. 



Table 5. Correlation coefficients for SA94 QSOs 
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Fig. 7. Variability index versus absolute magnitude for the 
QSOs in the SA94 (•), SGP (x) and SA57 (□) 
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Fig. 8. Variability index versus redshift for the QSOs in the 
SA94 (•), SGP (x) and SA57 (□) 



nitude and anticorrelated with redshift (correlation co- 
efficients respectively of +0.24 and —0.20, both with a 
significance greater than 99 per cent). More luminous 
and/or higher redshift objects are less variable, but, be- 
cause of the strong anti-correlation between absolute mag- 
nitude and redshift, it is difficult to ascertain if one of 
the two correlations is spurious. This result is substan- 
tially in agreement with previous works (see, for example, 
Cristiani et al. 1990). 

Then, we have merged the SA94 with the SA57 and 
SGP samples, and re-computed the correlation matrix 
(Table 6). 



Table 6. Correlation coefficients for the total of QSOs of the 
three samples 
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Mb 


IDX 


0.21 ± 0.04 


-0.12 ± 0.04 


0.27 ±0.04 


Mb 


0.39 ± 0.04 


-0.73 ± 0.02 
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0.29 ± 0.04 







The anticorrelation between absolute magnitude and 
redshift is reduced. The correlation of the variability index 
with Mb results unvaried and its significance is increased. 
The anticorrelation of the variability index with redshift 
is weaker but marginally significant. It appears now also a 
correlation with B (0.21 with significance > 99 per cent), 
a trend caused by the objects of the SA57 sample, less 
luminous and more variable, as we can see in Fig. 6. 

To further investigate the dependences of the variabil- 
ity index on the absolute magnitude and redshift, as in 
the previous sub-section, we have subdivided the L — z 
plane in a number of sub-areas and evaluated the average 
variability index in each of them. The result is shown in 
Fig. 9. 

Going from upper to lower panels within the various 
redshift slices, the anti-correlation between variability in- 
dex and absolute magnitude is again apparent. A weaker 
dependence is observed on the redshift. If we examine 
the average variability index for the objects with —25 > 
M B > -27, we find that below redshift 1 < IDX >= 
0.030 ± 0.009, while for z > 1 < IDX >= 0.069 ± 0.005, 
a difference significant at a 3.8(7 level. 

As in the analysis of the SF, the bin at < Mb >= 
—24.3 and < z >= 2.28, shows a much higher variability 
than expected from the general trend. Such a discrepant 
behaviour might be the indication of a bias in the selec- 
tion technique in favour of variable QSOs. Color selection 
techniques are based on the separation of the QSO can- 
didates from the stellar locus in a multi-color space (or 
plane) and a significant bias is expected to occur when 
(typically for 2.2 < z < 3.5) a QSO is located adjacent 
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Fig. 9. The average rest-frame variability index as a function 
of absolute magnitude and redshift. In the upper left corner 
of each panel the number of objects in the bin is shown. In 
the lower left the average absolute magnitude and in the lower 
right the average redshift are given. In the middle the average 
variability index and its la uncertainty is reported 



to, or coincident with, the volume occupied by the stellar 
locus: QSOs whose detection probabilities are very small, 
because of variability, may be scattered into regions where 
they are detectable (e.g. Warren et al. 1994). The bin at 
< Mb >= —24.3 and 2 < z < 3 is based on 4 objects per- 
taining to the SA57 that have been selected in a U — J , 
J — F plane (Koo &: Kron 1988) and are consequently ex- 
posed to the above-mentioned bias. A similar bias might 
also affect the objects in the SGP with 2.2 < z < 3.5, dis- 
covered with multi-colour methods (Warren et al. 1994). 
We have compared the variability indices of a subsample of 
the SGP with 2.2 < z < 2.7 and -28.0 < M B < -25.5 (40 
objects) with a corresponding subsample of the SA94 (18 
objects) that, being selected with objective prism data, is 
not affected by this bias. The average variability indices 
turn out to be 0.069±0.012 and 0.066±0.021, respectively, 
showing that for the SGP sample this bias is not playing 
any significant role. 

Another way of disentangling the dependence of the 
variability on the redshift from the anticorrelation with 
the absolute luminosity is given by the method of partial 
correlation analysis (Spiegel 1991). Applying this recipe to 
the results of Table 6, a value of 0.12 ±0.04 (99 per cent) 
results for the correlation coefficient between variability 
index and redshift. Again it appears reasonably well estab- 
lished that, for a given absolute luminosity, higher redshift 



(due to the redshift-absolute magnitude anticorrelation) 
at first glance of the Table 6. The results do not change 
if one of the samples is eliminated, they simply become 
less significant, showing that the different time sampling 
of the different samples does not lead to biases. 

The same analysis can be carried out in the observer's 
rest-frame, estimating the variability indices between 6 
and 15 yr. The correlation matrix turns out to be: 



Table 7. Correlation coefficients computed in the rest-frame 
of the observer 
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If we remove the dependence of the variability index on 
the absolute magnitude with the method of partial corre- 
lation analysis, the coefficient of correlation between vari- 
ability index and redshift becomes +0.02 ± 0.04. How- 
ever, if we examine the average variability index for the 
objects with —25 > Mb > —27, we find that below 
redshift 1 < IDX >= 0.043 ± 0.008, while for z > 1 
< IDX >= 0.086±0.008, a difference significant at a 
3.8(7 level. This latter result seems in contradiction with 
the analysis of the SF and of the correlation matrix. It 
has to be considered that it shows a peculiar behaviour 
of 17 objects at low redshift with respect to 234 objects 
at high redshift. The SF fitting or the correlation matrix 
analysis are of course insensitive to such a small subset in 
the L — z plane (in fact the low-z bins show a high % 2 per 
bin). 



5. Discussion and Conclusions 

By merging three different samples (SA94, SGP, SA57), 
we have been able to analyse the variability of a statisti- 
cally well defined set of 486 optically selected QSOs. 

The ensemble structure function and individual vari- 
ability indices of the QSOs have been evaluated by means 
of "robust" statistical estimators, less exposed to the influ- 
ence of possible "outliers" with respect to more conven- 
tional estimators. The effects of the photometric errors 
have been subtracted, allowing a meaningful comparison 
of the three different samples. 

Although the coverage of the Mb — z plane of the 
present samples is not completely uniform, the definition 
of suitable sub-samples has allowed us to disentangle the 
average relationships between variability and redshift and 
between variability and luminosity. 



structure function and the individual variability indices in 
the QSOs rest frame show that: 

1. A negative correlation between variability and lumi- 
nosity is clearly present, in the sense that more lu- 
minous QSOs show less variability (in magnitude), 
confirming previous results (Cristiani et al. 1990, 
Hook et al. 1994). 

2. A significant positive correlation exists between vari- 
ability and redshift, as suggested already by Giallongo 
et al. (1991). 

3. Such correlations may be equally well parameterized 
either with a model in which the timescale of the vari- 
ability is fixed for all the QSOs and the amplitude lin- 
early increases with the absolute magnitude and red- 
shift, or with a model in which the timescale of the 
variability linearly depends on the absolute magnitude 
and the amplitude is only a function of the redshift. 

The same analysis carried out in the observer's frame 
provides the following results: 

1. There is a negative correlation between variability and 
luminosity. 

2. The timescale of variability does not depend signifi- 
cantly neither on the absolute magnitude nor on the 
redshift. 

3. The ensemble structure function is well represented by 
a parameterization in which, with a fixed timescale of 
about 5.5 yr, the amplitude linearly increases with the 
absolute magnitude. 

4. Although the general behaviour of the SF does not 
show a systematic variation of the timescale and/or 
amplitude with redshift, if we examine the average 
variability index for objects with —25 > Mb > —27, 
we find that below redshift 1 quasars are significantly 
less variable than at higher redshift. 

The anti-correlation observed in the QSOs rest frame 
between absolute luminosity and variability is considered 
by some authors as an evidence in favor of the sub-units 
model, in which variability is caused by an ensemble of 
individual flares. To test this hypothesis, following Pica 
and Smith (1993), we have made a simple calculation. Let 
a QSO be composed by N random flaring sub-units; in this 
case the signal is proportional to N , whereas the noise (e.g. 
a supernova event in the starburst model) is proportional 
to VW. Then S/N k, N/VW = VN oc VI, where L is the 
luminosity of the QSO. Thus the amplitude of the relative 
variability should decrease as \ j\fTi and the SF as 1/L. 
A difference of 2 magnitudes, corresponding to a ratio 
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should give rise to an equal ratio between the amplitudes 
of the SF 
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value < 2, unacceptably lower than the predictions of the 
sub-units model, at least in its simplest form. However, 
any additional background luminosity contributing to the 
total luminosity of the QSOs would introduce departures 
from the L~ x relationship. This is, indeed, the case of 
the starburst model for AGNs (Terlevich et al. 1992), in 
which variability is caused by the supernova explosions 
of a nuclear young stellar cluster. The stellar background 
in this model accounts for about half the luminosity in 
B band (Aretxaga &: Terlevich 1994), and tends to flat- 
ten the L~ x law (Aretxaga 1993; Aretxaga et al. 1993). 
Even a more standard model of a black-hole which in- 
duces pulses of light (due to, for example, accretion events 
or stellar impacts onto the disk) would introduce some 
background luminosity in a quiescent stage that, poten- 
tially, could modify that law. Another important effect is 
expected to be due to the variability-redshift correlation. 
If such a correlation is not accounted for, since the more 
luminous QSOs are those of higher redshifts, the ampli- 
tude of their observed B-band variability would tend to be 
larger than what expected on the basis of the luminosity- 
variability anticorrelation, causing a flattening of the a 
priori L~ x law. This effect is important enough to in- 
troduce serious departures from the standard law for the 
sub-unit model, and account for the flattening observed 
(Cid Fernandes 1995). 

Comparisons with more refined models of variability 
produced by dense stellar clusters and supernova explo- 
sions are being developed, and will be discussed in detail 
in a separate paper. Unfortunately these models are the 
only ones that, so far, permit a detailed comparison be- 
tween observed light curves and the physical parameters 
of the events that produce the variability (energy, metal- 
licity, etc). 

The positive variability-redshift correlation may 
be interpreted as an increase of variability at the 
shorter and shorter wavelengths redshifted in the fixed 
observational bandpass (Giallongo et al. 1991), as al- 
ready observed in individual objects (Edelson et al. 1990, 
Kinney et al. 1991, Paltani &: Courvoisier 1994). Such a 
behaviour should be connected to the physical mecha- 
nisms governing the energy flux at different wavelengths. 
In order to confirm this result and its interpretation it 
will be important to study the ensemble variability of the 
same QSO samples on red plates. 

Hawkins (1993) has recently proposed that nearly all 
QSOs are being microlensed. According to him, since the 
redshift distribution of the lensing objects (for a uni- 
form population of lenses) is expected to strongly peak 
at z ~ 0.5, with only weak dependence on the redshift of 
the source, for such a model one would not expect to see 
a significant increase in time scales with redshift (in the 
rest-frame of the observer). The present data are not able 
to disprove such a behaviour. It remains to be established 
whether the absence of a dependence of the variability 



the structure function, is a true result of microlensing or 
the time dilation effect counterbalances the evolution of 
the intrinsic variability with redshift. On the other hand, 
it has also to be clarified if the significant decrease of the 
variability index of QSOs with —25 > Mb > —27 found 
at z < 1 is due to an intrinsic change of the QSO variabil- 
ity or to a smaller expected frequency of the microlensing. 
In this respect, it will be necessary to carry out further 
tests, for example on the achromaticity of the variability 
of individual objects, for which again the study of the en- 
semble variability of the same QSO samples on red plates 
will be extremely important. 
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A. Taking into account measurement errors in the 
computation of the structure function 

The difference between the magnitudes at two epochs, in pres- 
ence of errors, is given by 



m; + e; 



(Al) 



where m; is the true magnitude of an object at the i-th epoch 
and e;, ej are the errors. But, because of the triangular inequal- 
ity, 



■ m; I > |m; + e; — m; 



(A2) 



where the term at the first member is the difference of the "true 
magnitudes" we are looking for and cannot be calculated in a 
straightforward way. 

To subtract the influence of the errors we have adopted the 
following procedure: 

1. The various plates have been grouped in epochs accord- 
ing to the time of observation (12 for the SA94 sample, see 
Table 2; 11 for SA57; 7 for SGP). 

2. For each object, the median of all the plate-magnitudes 
has been calculated. To avoid biases, for each plate, the 
measurements of all the objects with a median magnitude 
fainter than the plate limit-magnitude have been excluded 
from the following computations. 

3. The epoch-magnitudes EMAG have been calculated as the 
mean of the magnitudes of each epoch. 



\EMAGj - EMAGj 



(A3) 



has been computed, where <r ; , aj are the magnitude un- 
certainties at the epochs i, j. The Uij have been assigned 
to the appropriate temporal bin, identified by 

Ai = 1**JlM (A4) 

when the SF is computed in the QSOs rest-frame or simply 
by 

Ai = \U - tj\ (A5) 

when the SF is computed in the observer'ss rest-frame. 
5. For each temporal bin the average deviation 



DEV(At) 



N 



(A6) 



has been evaluated, where the sum is carried out on the N 
magnitude differences pertaining to the bin. This quantity 
DEV(At) is the result of the magnitude variations due to 
the intrinsic QSO variability plus the photometric errors of 
each epoch-magnitude. 
6. We have then computed the quantity 



DEV s (At) 



N 



(A7) 



where \vij\ is the expected absolute deviation for the 
epochs i and j, assuming that the difference of the epoch- 
magnitudes of an object is the result of a gaussian process 
with a variance resulting from the sum of the variances 
due to the appropriate photometric uncertainties plus a 
variance SF e due to intrinsic variability. The correct value 
of the SF e in a given temporal bin is assumed to be the 
one for which 

DEV - DEV S = f(SF e ) = (A8) 

7. Since the differences \EMAGi — EMAGj \ are not indepen- 
dent measurements, to estimate the errors we have used 
the bootstrap sampling method (Barrow et al. 1984). The 
bootstrap method mimics the process of repeating the in- 
put samples a large number of times. To generate the artifi- 
cial samples we selected 486 times a QSO at random from 
the original samples. Thus the same QSO could appear 
more than once in the artificial sample. The process was 
repeated 100 times and for each artificial sample the SF e 
was computed. For each time bin of the SF e the RMS of 
the estimates obtained from the 100 simulations has been 
adopted as the uncertainty on the value computed in the 
Eq. (A8). 
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